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Abstract. 

The Modern-Era Retrospective Analysis for Research and Applications-2 (MERRA2) version of 
the GEOS-5 Atmospheric General Circulation Model (AGCM) is currently in use in the NASA 
Global Modeling and Assimilation Office (GMAO) at a wide range of resolutions for a variety of 
5 applications. Details of the changes in parameterizations subsequent to the version in the original 
MERRA reanalysis are presented here. Results of a series of atmosphere-only sensitivity studies are 
shown to demonstrate changes in simulated climate associated with specific changes in physical pa- 
rameterizations, and the impact of the newly implemented resolution-aware behavior on simulations 
at different resolutions is demonstrated. The GEOS-5 AGCM presented here is the model used as 
1 0 part of the GMAO’s MERRA2 reanalysis, the global mesoscale ’’nature run”, the real-time numer- 
ical weather prediction system, and for atmosphere-only, coupled ocean-atmosphere and coupled 
atmosphere-chemistry simulations. 

The seasonal mean climate of the MERRA2 version of the GEOS-5 AGCM represents a sub- 
stantial improvement over the simulated climate of the MERRA version at all resolutions and for 
1 5 all applications. Fundamental improvements in simulated climate are associated with the increased 
re-evaporation of frozen precipitation and cloud condensate, resulting in a wetter atmosphere. Im- 
provements in simulated climate are also shown to be attributable to changes in the background 
gravity wave drag, and to upgrades in the relationship between the ocean surface stress and the 
ocean roughness. The series of ’’resolution aware” parameters related to the moist physics were 
20 shown to result in improvements at higher resolutions, and result in AGCM simulations that exhibit 
seamless behavior across different resolutions and applications. 
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1 Introduction 


The various activities of NASA’s Global Modeling and Assimilation Office (GMAO) neccesitate a 
model that can function seamlessly across many different resolutions and applications. These appli- 
cations include real-time atmospheric analyses and forecasts at a resolution of 0.25°, long term re- 
analyses at 0.5°, coupled atmosphere-ocean and coupled atmosphere-chemistry simulations at 1-2°, 
and global mesoscale simulations at 7km. The Modern-Era Retrospective Analysis for Research 
and Applications-2 (MERRA2) version of the GEOS-5 AGCM is part of an ongoing development 
of a new generation AGCM at GMAO. The focus of the development of the MERRA version of 
the GEOS-5 AGCM was on the behavior of the components of the hydrological cycle in reanalysis 
mode, while the focus of the development of the MERRA2 AGCM was on a model that functions 
seamlessly in numerical weather prediction, reanalysis, climate and global mesoscale modes. To 
this end, some of the physical parameterizations were replaced, some parameters governing the 
behavior of other physical parameterizations were changed, and resolution-aware parameters were 
implemented in the moist process parameterizations. 

Many studies exist that describe major improvements in new versions of AGCMs and show the 
improvements in simulations as compared to reanalyses and other observations (ie., Neale et al., 
2013, Donner, et al., 2011, Pope et al., 2000). The present study adds to that type of analysis by 
carefully documenting the connection between individual changes in the physical parameterizations 
of the AGCM and improvements in the climate simulation at coarse resolution. A series of sen- 
sitivity experiments were conceived and analyzed to explore, step by step, each important change 
in parameterizations between the MERRA and MERRA2 AGCMs, and to demonstrate the impact 
on the simulated climate. The present study also describes and analyzes the improvements in high 
resolution simulations due to some changes in parameterizations specifically targeted for those res- 
olutions. 

The details of the changes in the AGCM physical parameterizations are described in the next 
section, the step by step experiments to isolate the impacts of these changes are described in section 
3, and the impacts of the “resolution aware” aspects of the AGCM are described in section 4. The 
study is summarized in section 5. 

2 Description of the MERRA2 version of the GEOS-5 AGCM 

The generation of the GEOS-5 Atmospheric General Circulation Model (GCM) that was used as 
part of NASA’s Modern-Era Retrospective Analysis for Research and Applications (MERRA) is 
described in Rienecker, et al. (2008), and most of the subsequent development of the physical 
parameterizations for the current, MERRA2 version is described in Molod et al. (2012). In addition 
to the changes in the physical parameterizations, the development of the MERRA2 AGCM also 
included two fundamental elements that will not be addressed in the present study. The horizontal 
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discretization of the MERRA2 AGCM is computed on the cubed sphere grid of Putman and Lin 
(2007), although it still retains the option to use the latitude/longitude discretization. The cubed 
sphere grid allows for the relative uniformity of grid spacing at all latitudes, and avoids the grid 
spacing singularities found in the latitude/longitude grid. In addition, the MERRA2 AGCM has been 
modified to account for the change in total mass due to the change in total water content computed 
in the moist and turbulence processes. The total mass of each layer is adjusted to include these 
changes in total water content, and the associated adjustment is then made to the specific masses 
of all constituents, including water substances. The benefit for the AGCM mean simulated climate 
is small, but results in the conservation of dry mass during the simulation. The algorithm for this 
adjustment and the benefits for AGCM simulations and data assimilation experiments are described 
in detail in Takacs et al. (2014). 

A brief summary of the model’s physical parameterizations relevant to the present study is pro- 
vided here. The GEOS-5 AGCM physics includes parameterization schemes for atmospheric con- 
vection, large scale precipitation and cloud cover, longwave and shortwave radiation, turbulence, 
gravity wave drag, a land surface model, a thermodynamic sea ice model, and a simple glacier 
model. 

Convection is parameterized using the Relaxed Arakawa-Schubert (RAS) scheme of Moorthi and 
Suarez (1992) and includes a scheme for the generation and re-evaporation of falling rain (Bacmeis- 
ter et al., 2006). A ’’stochastic Tokioka trigger” function (Bacmeister and Stephens, 2011) governs 
the upper limits on the allowable entrainment by sampling from a probability distribution function 
with specified parameters. The prognostic cloud cover and cloud water and ice scheme is from 
Bacmeister et al. (2006). The scheme includes large scale condensation governed by the probability 
distribution function described in Molod (2012), evaporation, autoconversion and accretion of cloud 
water and ice, sedimentation of cloud ice and re-evaporation of falling precipitation. 

The turbulence parameterization is based on the non-local scheme of Lock (2000) scheme, acting 
together with the Richardson-number based scheme of Louis and Geleyn (1982). The original Lock 
scheme was extended in GEOS-5 to include moist heating and entrainment in the unstable surface 
parcel calculations. The Monin-Obukhov similarity theory based parameterization of surface layer 
turbulence is described in Helfand and Schubert (1995), and includes the effects of a viscous sublayer 
for heat and moisture transport over all surfaces except land. The ocean roughness is determined 
by a polynomial which is a blend of the algorithms of Large and Pond (1981) and Kondo (1975), 
modified in the mid-range wind regime based on recent observations in the southern ocean according 
to Garfinkel et al. (201 1) and in the high wind regime according to Molod et al. (2013). 

The longwave radiative processes are described by Chou and Suarez (1994), and the shortwave 
radiative processes are from Chou and Suarez (1999). The gravity wave parameterization computes 
the momentum and heat deposition into the grid-scale flow due to orographic (McEarlane, 1987) and 
nonorographic (after Garcia and Boville, 1994) gravity wave breaking. The Land Surface Model 
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from Koster et. al (2000) is a catchment-based scheme that treats subgrid scale heterogeneity in 
surface moisture statistically. Glacial thermodynamic process are parameterized using an adaptation 
of the Stieglitz et al. (2001) snow model to glacial ice (Cullather et ah, 2014), and the catchment 
and glacier models are each coupled to the multi-layer snow model of Stieglitz et al. (2001). Sea ice 
albedos in the northern hemisphere are from the monthly mean observations of Duynkerke and de 
Roode (2001). 

3 Evolution of Low Resolution Simulated Climate from MERRA AGCM to MERRA2 AGCM 

The mean climate characteristics of a single 30-year MERRA2 AGCM simulation on the lati- 
tude/longitude grid at a horizontal resolution of 0.5° x 0.5° were evaluated by comparison with 
reanalysis and with different satellite and in situ based observational estimates (Molod et ah, 2012). 
They found substantial improvements in some key aspects of the mean circulation in the MERRA2 
version of the GEOS-5 AGCM, and also reported on existing discrepancies between the modeled and 
observed climates. Here we present the results of a series of experiments designed to attribute each 
fundamental improvement in AGCM simulated climate to a specific change in parameterization. 

The experiments to be described in this section were all conducted on the latitude/longitude grid 
at 2° X 2.5° horizontal resolution, on a vertical hybrid eta-pressure coordinate grid with 72 levels, 
spaced to increase the resolution near the surface and near the tropopause. The simulations were 
all forced with observed sea surface temperatures (Reynolds, 2002), and ran for 30 years each. 
The sequence of experiments was designed to start with the MERRA2 AGCM as the control and 
backtrack, one parameterization change at a time or small groups of parameterization changes as 
a time, to a model that replicates the MERRA AGCM simulated climate. The parameterization 
changes are listed in table 1, and the full sequence of the control and 7 sensitivity experiments to be 
described in this section is listed in Table 2. 

3.1 Ocean Surface Winds 

The parameterization of the surface layer turbulence in the MERRA2 AGCM includes a substantial 
modification of the functional relationship between ocean surface roughness and wind stress, shown 
in Eigure 1 . The relationship for the moderate range surface wind speeds from the MERRA AGCM 
(green) and the MERRA2 AGCM (black) is shown in figure la, where the increased roughness 
based on the implementation of Garfinkel et al. (2011) is apparent. The first experiment in the 
series, experiment 1, reverts back to the formulation for the relationship between ocean surface 
roughness and stress used in the MERRA AGCM. The effect for simulations at 2.° x 2.5° resolution 
is expected be an increase in the experiment 1 simulated surface wind speeds in the mid-range of 
wind speeds, that is, in the 5 m to 25 m range. Figure lb shows the relationship between 
wind speed and roughness for a larger range of wind speeds after Molod et al. (2013), where the 
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reduction of roughness at speeds greater than approximately 30 m is apparant. This change is 
expected to result in a net increase of wind speeds in higher wind regimes as was shown in Molod et 
al. (2013) for simulations at 0.25 ° resolution, but this impact is not apparent at the resolution of the 
experiments described here because the simulated wind speeds generally do not reach 30 m 

Surface wind speeds from the MERRA2 AGCM control and experiment 1 are shown in figure 2. 
The change in the simulated surface winds is most apparent in the southern hemisphere, where the 
Goddard Satellite-based Surface Turbulent Fluxes (GSSTF, Shie et ah, 2009) surface winds, seen in 
figure 2b, show values near 8 m s“^, experiment 1 (2d) shows surface winds near 12 m s“^ and the 
MFRRA2 AGCM control (2a) shows surface winds near 10 m s“^. The difference from the GSSTF 
estimate (shown in 2c and 2f) shows a reduction in the difference from up to 4 m s“^ in experiment 
1 to a difference of up to 2 m s“^ in the MFRRA2 AGCM control, pointing out the improvement in 
AGCM simulated climate due to the change in roughess formulation. 

3.2 Quasi-biennial Oscillation 

The latitudinal profile of background nonorographic drag in the MFRRA2 AGCM was modified to 
include a specified source related to tropical precipitation in addition to the local maxima related to 
storm track precipitation. The background profiles used in the MFRRA2 and the MFRRA AGCMs 
are shown in figure 3. Experiment 2 of the series examines the impact of the change in the gravity 
wave drag parameterization of background drag in the tropics, and returns to the MFRRA AGCM 
background drag. Experiment 2 is therefore expected to exhibit a stratospheric wind with no quasi- 
biennial oscillation (QBO) variability. The tropical (latitude range lOS to ION) zonal average of the 
zonal wind as a function of height and time is shown in figure 4. The patterns of large easterly and 
westerly winds that slant downwards in pressure as time proceeds indicate the downward propagation 
of the variations due to the QBO. The QBO pattern is seen in figure 4a and 4c, which are the results 
from experiment 1 from and MFRRA reanalysis (Rienecker, et ah, 2011), respectively. Figure 4b, 
however, that shows the results of experiment 2 that uses the old background drag formulation (from 
figure 3), shows no QBO pattern of variability. 

3.3 Stable Surface Layer Fluxes 

The parameterization of the surface layer in the MERRA2 AGCM uses the scheme of Helfand and 
Schubert (1995) based on Monin-Obhukov (MO) similarity theory. The scheme replaced the Louis 
(1979) scheme used in the MFRRA AGCM. The implementation of the MO scheme included the 
use of a different set of stable layer stability functions, and a different formulation for the viscous 
sublayer (the laminar layer that can act to impede the flux of heat and moisture). The stable surface 
layer stability functions in the Helfand scheme result in an increased turbulent heat exchange (of 
both signs) under stable conditions. Figure 5 shows a scatter diagram of the sensible heat flux 
as a function of surface bulk Richardson number under conditions where the monthly mean air 
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temperature exceeds the monthly mean skin temperature. The black points are from the simulation 
with the Helfand surface layer, and the red are from the simulation using the Louis scheme. The 
larger values of sensible heat flux in the Helfand simulation are apparent, and even more apparent 
when the monthly mean sensible heat flux is downward. 

Experiment 3 was designed to examine the impact of the change in the surface layer parame- 
terization by reverting back from the Helfand and Schubert scheme to the Louis scheme. Figure 
6 shows the sensible heat flux from experiments 2 and 3 along with the difference between them. 
The differences shown here are attributable to the removal of the viscous sublayer over land in the 
Helfand and Schubert scheme and to the change in the stable layer stability functions. Figure 6c 
shows that over most land surfaces the difference in sensible heat flux is negative, indicating less 
sensible heating when using the Louis scheme. This sign of the difference is consistent with the 
removal of the viscous sublayer over land surfaces in the Helfand scheme, that would remove some 
resistance to turbulent exchange that is present over unvegetated land surfaces in the Louis scheme 
due to the viscous sublayer. Figure 6c also shows regions where the sensible heat flux is greater in 
the Louis scheme than in the Helfand scheme. These are regions where the surface layer is stable, 
and where the sensible heat itself is largely downward (that is, the air temperature is greater than the 
skin temperature). The change in stability functions between the Louis and Helfand scheme, that 
allows more turbulent exchange in the Helfand scheme, is consistent with the sign of the difference 
in sensible heat flux in regions where the heat flux is downward. 

3.4 Atmospheric Moisture, Clouds and Stationary Wave Pattern 

3.4.1 Critical Relative Humidity 

The algorithm for large scale condensation, as described in Bacmeister et al. (2006), assumes that 
the probability distribution function (PDF) of total water is “top-hat” shaped. The width of the 
PDF can be shown to be associated with a “critical relative humidity” (RHcrit) that governs cloud 
macrophysical and microphysical processes such as condensation and evaporation (Molod, 2012). 
The relationship between RHcrit and PDF width is such that a wider PDF corresponds to a lower 
RHcrit- The MERRA2 AGCM RHcrit (Molod, 2012) represents a change in both the magnitude 
and vertical structure from the RHcrit in the MERRA AGCM. Typical RHcrit profiles from the 
MERRA and MERRA2 AGCMs are shown in figure 7, and indicate generally lower values in the 
MERRA2 AGCM formulation except in the boundary layer, where turbulent mixing is sufficient to 
homogenize the total water distribution and so result in narrower PDF. 

Experiment 5 was designed to examine the impact of the change in RHcrit- Removing this 
change, which for much of the atmosphere means a larger RHcrit, should result in a simulation 
that is generally wetter because the atmosphere is being adjusted back to a higher relative humidity 
(RH). The zonal mean relative humidity from experiments 4 and 5, along with the difference between 
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them, is shown in figure 8. The experiment 5 minus experiment 4 difference shows a clear increase 
in relative humidity in the MERRA AGCM-like experiment due to the increase in RHcrit ■ Relative 
to available observational verification, the MERRA2 AGCM shows a general wet bias (Molod et 
al., 2012), which means that the experiment 4 relative humidity field is closer to the observed than 
experiment 5’s RH field. 

In addition to having a substantial impact on atmospheric moisture, the change in RHcrit also had 
an impact on the distribution of cloud cover. The higher RHcrit simulation (the MERRA AGCM- 
like experiment) could either be expected to have less cloud cover because the atmosphere must 
contain more moisture before new cloud water will be condensed, or could be expected to have more 
cloud cover due to the feedback of a generally moister atmosphere. Figure 9 shows the zonal mean 
cloud cover from experiments 4 and 5, and an observational estimate of zonal mean cloud cover 
from AIRS. The MERRA AGCM-like experiment, experiment 5 (figure 9a), shows increased cloud 
cover in the 300-600 mb range relative to the MERRA2 AGCM-like experiment, in particular at 
high latitudes in both hemispheres. In this regard, the MERRA2 AGCM-like experiment result more 
closely resembles the AIRS cloud cover estimate (figure 9c). This result confirms the hypothesis 
that higher RHcrit results in an atmosphere that is wetter in the mean and has more cloud. The 
MERRA-like AGCM result (experiment 5) also shows smaller cloud cover near the boundary layer 
at almost all latitudes. At high latitudes, the MERRA2 AGCM-like boundary layer cloud more 
closely resembles the AIRS estimate, in the tropics the MERRA2 AGCM-like boundary layer cloud 
is larger than AIRS, while the MERRA AGCM-like boundary layer cloud is smaller. The change in 
boundary layer cloud between the MERRA AGCM-like experiment and the MERRA2 AGCM-like 
experiment is not consistent with the free atmosphere response to the RHcrit change because the 
presence of boundary layer turbulence makes RHcrit less of a determining factor for model mean 
relative humidity there. 

3.4.2 Re-evaporation of Precipitation and Condensate 

The MERRA2 AGCM scheme for the re-evaporation of precipitation and suspended cloud water and 
ice contains a series of new parameter settings that result in a substantial increase over the MERRA 
model in the re-evaporation of snow and ice. The impact of the changes in parameter settings on the 
water vapor source due to re-evaporation for the December- January-February average is shown in 
figure 10. The largest increase in the MERRA2 AGCM is aloft, near 500 mb, where the increase is 
up to 0.7 g kg~^ day~^. 

Experiment 6 examines the impact of the change in re-evaporation of snow and ice in the MERRA2 
AGCM, that is perhaps the most crucial parameterization change. The removal of this change is ex- 
pected to result in a drier atmosphere, in particular aloft. Because of the importance of this change in 
parameterization, this simulation is expected to resemble in large part the climatology of the AGCM 
used as part of MERRA. Figure 1 1 shows the direct impact of the change from the MERRA2 AGCM 
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to the MERRA AGCM re-evaporation, and, as expected, shows the drying related to the reduced re- 
evaporation in the MERRA AGCM-like experiment. Figures 1 la-c show the change in total precip- 
itable water, and the difference (MERRA AGCM formulation minus MERRA2 AGCM formulation, 
shown in 1 Ic) is always negative everywhere. The vertical distribution of the moisture is shown with 
the specific humidity in figures lld-f, where the difference plot (Ilf) also shows an almost global 
reduction in atmospheric water vapor. 

The resulting mean circulation in boreal winter underwent a substantial change associated with 
this drying, and represents the most substantial impact on the simulated climate of all the elements 
of the MERRA AGCM to MERRA2 AGCM transition. Figure 12 shows the subtantial impact 
that changing the moisture levels had on the 30-year averaged December- January-February total 
precipitation. Given the Global Precipitation Climatology Project (GPCP, Huffman et al., 1995) 
climatology as a reference, the MERRA2 AGCM-like experiment (experiment 5, figure 12d) exhibits 
an Intertropical Convergence Zone (ITCZ) structure that is properly placed in longitude, a South 
Pacific Convergence Zone (SPCZ) that properly slants from the western Pacific to the southeast, and 
more accurately simulated storm track maxima. The precipitation field represents improvements 
in all these areas related to the MERRA AGCM-like simulation (experiment 6) shown in figure 
12a. The change in the tropical precipitation field related to the additional drying in experiment 6 
also resulted in substantial changes in the Pacific teleconnection patterns, specifically the Pacific- 
North- America (PNA) pattern in boreal winter. The eddy height field is an indicator of the strength 
of the PNA, and is shown in figure 13 in relation to the eddy height from MERRA reanalysis. 
The MERRA2 AGCM-like simulated PNA pattern (figure 12d) has a stronger and more properly 
oriented ridge near the west coast of North America relative to the PNA as simulated by the MERRA 
AGCM-like experiment (figure 12a). This change has implications for the poleward propagation of 
heat and momentum. The standard deviation of the difference from MERRA reanalysis estimates 
is also substantially reduced in experiment 5 (16.8 m) relative to experiment 6 (20.3 m). The direct 
impact of the change in re-evaporation was also evident in the boreal summer climatology. Figure 14 
shows this both in the total precipitable water fields and in the specific humidity fields, that exhibit 
differences between experiments 5 and 6 that are of the order of the differences seen in boreal winter. 
The impact on the mean summertime circulation, however, was minimal. 

3.5 Breakup of the Southern Hemisphere Stratospheric Jet 

The changes to the gravity wave drag parameterization included the modification of the ’’intermit- 
tency factor”, used to reduce the strength of the gravity wave drag based on expected departure from 
linear theory. The value of the intermittency factor was increased in the MERRA2 AGCM for oro- 
graphic waves as a function of latitude, changing from a MERRA AGCM global value of 0.125 to 
values reaching 0.3125 south of approximately 40°S. 

Experiment 7 is the last experiment in the series, incorporating the effects of all the fundamental 
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changes between the MERRA and MERRA2 versions of the GEOS-5 AGCM. The removal of the 
intermittency change is expected to decrease the orographically induced drag in the southern hemi- 
sphere, thereby depositing less momentum aloft and increasing the strength of the westerlies. The 
zero-wind contour in the southern hemisphere can be used as an indicator of the level and timing of 
the stratospheric jet breakup. Figure 15 shows the zero-wind contour from experiment 7, experiment 
6, and MERRA. At the 2 degree horizontal resolution of the simulations described here, the zero 
wind line is higher in altitude and delayed in time relative to MERRA in both simulations, reflecting 
a delayed stratospheric jet breakup. The decreased intermittency factor in experiment 7, however, 
delays the jet breakup even more, showing the improvement in the MERRA2 AGCM-like simulation 
relative to the MERRA AGCM-like simulation. At higher spatial resolution (not shown) the increase 
of intermittency factor in the MERRA2 AGCM is effective in producing a reasonable evolution of 
the polar vortex breakdown in the southern Hemisphere. 

4 High resolution simulations and resolution aware behavior in the MERRA2 AGCM 

The modifications to the MERRA2 AGCM physical parameterizations described above resulted in 
improvements in simulated climate at all the resolutions relevant to GMAO. Additional develop- 
ments were implemented in the MERRA2 AGCM that were particularly applicable to higher res- 
olution (0.25° or higher) simulations. These included the changes in ocean surface roughness at 
high wind speeds (mentioned in section 3.1 and examined in detail in Molod et ah, 2013) and the 
implementation of ’’resolution aware” parameters. 

The implementation of the “resolution aware” behavior of the moist processes in the MERRA2 
version of the GEOS-5 AGCM was designed to improve the behavior of the high resolution simula- 
tions and to ensure more uniformity of model mean state across resolutions and applications. Figure 
16 shows an example of the undesireable behavior in the MERRA version of the AGCM that the 
’’resolution aware” parameters were meant to address. Panels 16a-c show the specific humidity at 
l.°, 0.5° and 0.25° resolution, and panels 16g-i show the increase of the bias in atmospheric mois- 
ture content relative to EC-Interim reanalysis (Dee et ah, 201 1) as the resolution increases. In order 
to mitigate this change in simulated climate with increasing resolution, the MERRA2 AGCM moist 
physics was modified to include two governing parameters that are specified a priori as a function 
of horizontal resolution. These are the critical relative humidity used for large scale condensation, 
and a parameter which governs the minimum allowable entrainment used for the ’’stochastic Tokioka 
trigger” of the convective parameterization. 

The dependance of RHcrit on horizontal resolution is based on the Molod (2012) analysis of a 
global mesoscale simulation, and is such that the RHcrit increases with finer resolution, as was seen 
in the progression from the 2° resolution curve up to the 0.25° curve in figure 7. This progression 
is consistent with an intuitive expectation that the variability of total water within an AGCM grid 
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cell decreases as the grid cell becomes smaller. The implementation of the horizontal resolution 
dependance of RHcrit (as described in section 2) in the MERRA2 AGCM resulted in an atmospheric 
moisture field that is more consistent across different resolutions, as seen in figure 17. 

The MERRA2 version of the GEOS-5 AGCM also includes a horizontal resolution dependant 
and stochastic Tokioka (1988) type trigger (described in Bacmeister and Stephens, 201 1) as part of 
the RAS convective parameterization. The cloud model in RAS computes the effect of individual 
entraining cloud plumes, and the trigger acts to effectively eliminate any cloud plume with too small 
an entrainment during ascent. Bacmeister and Stephens (2011) examined the observed relationship 
between the neutral bouyancy level of a particular sounding and the observed condensate (a proxy for 
the convective detrainment level), and found that this observational relationship can be approximated 
by sampling the minimum entrainment from a power law PDF. The parameters of the PDF are 
specified a priori, and vary with resolution. The stochasticity is designed to only occasionally permit 
the least entraining and therefore the deepest detraining cloud plumes. Lim et al., (2014) reported 
on a series of sensitivity studies to examine the impact of different choices of the PDF governing 
parameter on the simulation of strong tropical storms in the GEOS-5 AGCM. 

The resolution dependance of the stochastic Tokioka trigger is such that it more severly limits the 
parameterized convective mass flux at high spatial resolution, where we expect the larger scale con- 
vective updrafts to be resolved explicitly, and has little impact at low resolution. The change with 
resolution of the PDF parameters is shown in figure 18, presently specified in an ad hoc manner. 
Figure 19 shows a sequence of the resulting June averaged convective mass fluxes from simulations 
with different horizontal resolutions. The decrease of parameterized mass flux with resolution re- 
flects the increasingly restrictive trigger, selected from the PDFs with increasingly higher mimimum 
entrainment values shown in figure 18. The effect of this repressed RAS mass flux on the simulated 
climate is reflected in the total change in moisture due to moist processes, shown in figure 20. The 
cloud model in RAS includes a grid scale subsidence to compensate for the updraft mass flux that 
results in a drying of the sub-cloud layer. Figure 20 shows the decrease of the low level drying 
with increased horizontal resolution. The benefits of the reduced low level drying and the related 
maintenance of the cumulus available potential energy during tropical storm development at high 
resolution was demonstrated in the study of Lim et al. (2014). 

5 Conclusions 

Synthesis and Discussion 

The version of the GEOS-5 AGCM used in the Goddard Modeling and Assimilation Office 
(GMAO) MERRA2 reanalysis was developed for use across many different resolutions and applica- 
tions. A unique series of AGCM simulations were performed with the GEOS-5 AGCM to detail the 
impact of each change in parameterization between the MERRA version and the current MERRA2 
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version. The series of sensitivity experiments began with the current AGCM version and regressed, 
one step of development at a time, to the MERRA AGCM. 

The most substantial positive impact on the simulated AGCM climate was shown to be attributable 
to the increase of the re-evaporation of frozen cloud water and precipitation in the MERRA2 AGCM. 

The resulting atmosphere had a higher moisture content, and many aspects of the boreal winter 
climate were substantially improved relative to reanalysis. The moisture and cloud cover amounts 
were shown to be further improved by the implementation of an AIRS-based PDF of total water. 

The development of the MERRA2 AGCM also included the implementation of a set of parame- 
ters governing moist processes that contain an a priori change in behavior with horizontal resolution. 

The parameters are ones which govern the mimimum allowable entrainment into a convective up- 
draft, and govern the atmospheric relative humidity needed for the onset of condensation. Results 
of a limited set of experiments were shown to demonstrate the benefits of this “resolution aware” 
behavior at higher resolution. 

This study was focused on the results of atmosphere only simulations, but the resulting model has 
also performed well in coupled atmosphere ocean, coupled atmosphere chemistry, data assimilation, 
numerical weather prediction and global mesoscale applications. 

6 Code Availability 

The GEOS-5 source code is available under the NASA Open Source Agreement at http://opensource.gsfc.nasa.gov/projects/GEOS- 
5/. 
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Fig. 1. Scatter diagram of surface wind speed (msec versus ocean roughness (m) in the MERRA and 
MERRA2 AGCMs. a) diagram focusing on medium range of wind speeds, MERRA2 shown in black and 
MERRA in green, and b) diagram extending to high wind regimes, MERRA shown in green and MERRA2 in 
black. 
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Fig. 2. 30-year average December-January-February surface wind speed (m sec from: a) MERRA2 control, 
b) GSSTF, c) Control-GSSTF, d) Experiment 1, e) GSSTE, f) Experiment 1-GSSTF 
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Fig. 3. Background nonorographic drag from the MERRA and MERRA2 AGCM simulations 
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Fig. 4. Spatial average of zonal wind in m s ^ from lOS to ION latitude as a function of pressure level in mb 
and time from a) MERRA2 control, b) Experiment 2 and c) MERRA reanalysis. 
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Fig. 5. Surface bulk Richardson number as a function of sensible heat flux in W m 
experiment 3 (red) and experiment 2 (black). 


^ in a single July from 
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Fig. 6. June-July-August averaged sensible heat flux in W m ^ from a) Experiment 3 (Louis scheme), b) 
experiment 2 (Helfand scheme) and c) the difference, experiment 3 minus experiment 2. 
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Fig. 7. Critical relative humidity. Black from MERRA AGCM formulation, green from MERRA2 AGCM 
formulation for 1° resolution, and red from MERRA2 AGCM formulation for 2° resolution. 
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Fig. 8. December-January-February averaged relative humidity in percent from a) Experiment 5 (MERRA 
AGCM-like), b) experiment 4 (MERRA2 AGCM-like) and c) the difference, experiment 5 minus experiment 4. 
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Fig. 9. December-January-February averaged cloud fraction from a) Experiment 5 (MERRA AGCM-like), b) 
experiment 4 (MERRA2 AGCM-like) and c) AIRS retreivals 
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Fig. 10. The difference (MERRA2 AGCM-like minus MERRA AGCM-like) of zonal mean specific humidity 
source term due to all re-evaporation for December- January-Eebruary. 
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Fig. 11. a) December- January-Febmary total precipitable water in mm from experiment 6 (MERRA AGCM- 
like), b) same as a) but from experiment 5, c) same as a) but the difference, experiment 6 minus experiment 
5, d) December-January-February specific humidity in gkg“^ from experiment 6 (MERRA AGCM-like), e) 
same as d) but from experiment 5, f) same as d) but the difference, experiment 6 minus experiment 5. 
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Fig. 12. December- January-Febmary total precipitation in mm day ^ from a) experiment 6 (MERRA AGCM- 
like), b) GPCP, c) the difference, experiment 6 minus GPCP, d) experiment 5 (MERRA2 AGCM-like), e) GPCP, 
and f) the difference, experiment 5 minus GPCP. 
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Fig. 13. December- January-Febmary 300 mb Eddy Height Climatology in m from a) experiment 6 (MERRA 
AGCM-like), b) MERRA, c) the difference, experiment 6 minus MERRA, d) experiment 5 (MERRA2 AGCM- 
like), e) MERRA, and f) the difference, experiment 5 minus MERRA. 
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Fig. 14. a) June-July-August total precipitable water in mm from experiment 6 (MERRA AGCM-like), b) 
same as a) but from experiment 5, c) same as a) but the difference, experiment 6 minus experiment 5, d) June- 
July-August specific humidity in g kg“^ from experiment 6 (MERRA AGCM-like), e) same as d) but from 
experiment 5, f) same as d) but the difference, experiment 6 minus experiment 5. 
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Fig. 15. 30-year averaged annual cycle of zonal mean zonal wind, averaged from 70S to 50S latitude in m sec 
from Experiment 7 in red (MERRA AGCM-like), experiment 6 in blue (MERRA2 AGCM-like) and MERRA 
in green 
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Fig. 16. 
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Fig. 17. 
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Fig. 18. The probability distribution function for the mimimum entrainment allowed by the cumulus parame- 
terization for different AGCM horizontal resolutions. Black line is for 2 degree, blue for 1 degree, red for 1/2 
degree and green for 1/4 degree. 
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Fig. 19. 
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Fig. 20. 
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Table 1. Changes in GCM Algorithms from MERRA to Fortuna 


Module 

Algorithm Change 

Comments 

Moist 

Increased re-evaporation of precipitation 

Fundamental change in model climate 


Modified autoconversion 

Fundamental change in model climate 


Modified effective radius of cloud drops 

Fundamental change cloud forcing 


Anvil fractions cut in half 

Fundamental change cloud forcing 


Autoconvert ‘warm fog’ 

Important change in coupled simulations 


New critical RH with resolution dependance 

Substantial change in simulated moisture 


Cloud base set at PEL depth 

Remove clouds detraining below PEL height 


RAS time scale no longer depends on turbulence 



Stochastic RAS with resolution dependance 

Substantial impact at high resolution 


Turb 

Remove restrictions on diffusion from Louis 

Increase near surface diffusion 


Reformulate turbulent length scale in Louis 



Reduce Lock scheme when there is wind shear 

Impact on marine PEL 


Reduce cloud top entrainment for Lock plumes 

Impact on marine PEL 


Surf 

Implement Helfand and Schubert scheme 



Remove viscous sublayer over land surfaces 

Improve land temperatures 


Change ocean roughness for middle wind regimes 

Reduce wind bias in S. Ocean 


Change ocean roughness high wind regimes 

Increase tropical cyclone intensity 


Land Surf 

Changed parameters for evapotranspiration 

Impact on ratio of surface to canopy evaporation 


GW Drag 

Changed profile of background drag 

Substantial impact on QEO 


Added intermittency of drag 

Impact on timing of winter jet breakup 
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Table 2. Experiments to attribute MERRA to Fortuna AGCM simulation changes to changes in parameteriza- 
tions. 


Experiment 

Description 


Control 


Exp 1 

Back off change in ocean roughness 

Exp 2 

Exp 1 + Back off gravity wave background drag and surface hydrology 

Exp 3 

Exp 2 + use old surface layer parameterization 

Exp 4 

Exp 3 + Back off increase of Richardson-number diffusion 

Exp 5 

Exp 4 + Back off decrease of critical RH aloft and decrease below 

Exp 6 

Exp 5 + Back off increase of all re-evaporation 

Exp 7 

Exp 6 + Back off gravity wave drag intermittency 
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